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^  ABSTRACT  (jLJi 

The  report  presents  p  low-Froude-number  asymptotic 
expansion  for  the  far-field  wave-amplitude  function  defined  within 

the  Neumann-Kelvin  theoretical  framework#  This  diymptetie^ _ 

expansion  provides  a  simple  analytical  approximation  defined 
explicitly  in  terms  of  the  geometrical  characteristics  of  the  ship 
hull  and  the  disturbance  velocity  vector,  ffitp  low-Froude-number 
analysis  ffresen  ted  in-tbis-reporp shows  that  the  wave  resistance 
and  (hey far-field  wave  pattern  of  a  shipfstronzlv  [depend bn  the 
shape  of  the  hull,  notably  the  presence  of  flare  and  the  shape  of 
the  waterline  at  the  bow  and  tk$»stern.  fn  particu latythfy analysis 
predicts  that  the  nondimensional  wave-resistance  coefficient  is 
O(F0),  where  F  is  the  Froude  number,  for  a  ship  form  with  a 
y  region  of  flare,  0(F®)  for  a  ship  form  that  is  wall  sided 
everywhere  and  has  either  a  bow  or  a  stern  (or  bothj  that  is 
neither  cusped  nor  round,  and  0(F&)  for  a  wall-sided  ship  form 
with  both  bow  and  stern  that  are  either  cusped  or  round. 
^-The^qpalysis  also  shows  that  the  relative  importance  of  the 
nonlinear  terms  in  the  free-surface  boundary  condition  depends 
on  the  shape  of  the  hull.  Specifically,  the  contribution  of  the 
nonlinear  terms  in  the  free-surface  boundary  condition  to  the  far- 
field  wave-amplitude/unction,  K,  is  found  to  be  O(I0j[ 
irrespectively  of  fh^hull  form.jjhis  contribution  must  be 
compared  with  the  result  predicted  by  the  usual  Neumann-Kelvin 
theory,  in  which  the  free-surface  boundary  condition  is  linearized. 
This  linear  theory  predicts  that  K  is  0(F)  for  a  ship  form  having 
a  region  of  flare,  0(F2)  for  a  ship  form  that  is  wall  sided 
everywhere  and  has  either  a  bow  or  a  stern  (or  both)  that  is 
neither  cusped  nor  round,  and  0(F3)  for  a  wall-sided  ship  form 
with  both  bow  and  stern  that  are  either  cusped  or  round.  The 
contribution  of  the  nonlinear  terms  in  the  free-surface  boundary 
condition  is  then  asymptotically  negligible  for  ship  forms  having  a 
region  of  flare  or  a  sharp  wedge-like  bow  or/and  stern,  as  is  the 
case  for  a  majority  of  ships  (including  naval  ships). 
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ADMINISTRATIVE  INFORMATION 

This  study  was  funded  under  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center’s  Independent  Research  Program,  Program  Element  61152N, 
Project  Number  ZR02301,  Thsk  Area  ZR0230101,  Work  Unit  1542-108. 

INTRODUCTION 

The  wave  resistance  of  a  ship  advancing  at  constant  speed  in  calm  water  and  the 
amplitude  of  the  transverse  and  divergent  waves  in  its  far-field  wave  pattern  are 
defined  in  terms  of  the  far-field  wave-amplitude  function,  K(t),  by  means  of  simple 
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expressions,  which  may  be  found  in  Noblesse  [1,2]  for  instance.  The  function  K(t)  is 
given  by  the  sum  of  two  integrals,  namely  a  surface  integral  over  the  ship’s  mean 
wetted-hull  surface  and  a  line  integral  along  its  waterline,  involving  the  disturbance 
velocity  potential  in  their  integrands.  An  approximate  expression  for  the  function  K(t), 
valid  for  sufficiently  small  values  of  the  Froude  number,  was  recently  obtained  in 
Noblesse  [3]  by  approximating  the  hull-surface  integral  by  a  line  integral  along  the 
ship  waterline  and  combining  the  latter  waterline  integral  with  that  already  present  in 
the  original  unapproximated  expression  for  the  function  K(t).  The  resulting  low- 
Froude-number  approximate  expression  in  terms  of  a  waterline  integral  obtained  in  [3] 
is  now  given. 

Nondimensional  coordinates  (x,y,z)  =  (X,Y,Z)/L,  where  L  is  the  length  of  the 
ship,  attached  to  the  moving  ship  are  defined.  The  z  axis  is  taken  vertical  and 
pointing  upwards  and  the  x  axis  is  chosen  in  the  ship  centerplane  (port-and-starboard 
symmetry  is  assumed  here)  and  pointing  towards  the  bow.  The  flow  disturbance_due 
to  the  ship  is  represented  in  terms  of  the  nondimensional  velocity  potential  ^  =  <t>/UL, 


where  U  is  the  speed  of  the  ship. 

The  positive  half  of  the  mean  waterline  is  represented  by  the  parametric 
equations 

x  =  Xq(A)  and  y  =  y0(A) ,  (la,b) 

where  the  parameter  A  varies  between  its  bow  and  stern  values,  that  is 

Ag  ^  A  ^  A^  .  (1c) 

In  the  vicinity  of  the  mean  free  surface,  the  positive  half  of  the  hull  surface  is 
represented  by  the  parametric  equations 

x  =  Xq( A)  +  zxj(A)  +  z2x2(A)  +  ...  ,  (2a) 

y  =  yo(*)  +  zyi(*)  +  z2y2(*)  +  - .  (2b) 

where  AB  <  A  <  As  and  z  <  0  .  (2c) 

The  velocity  potential  <KA,z)  on  the  hull  surface  in  the  vicinity  of  the  plane  z  =  0 
likewise  is  expressed  in  the  form 

4  =  $o(A)  +  z+i(A)  +  +  —  •  (3) 

Differentiation  of  the  functions  xn(A),  yn(A),  4n(A)  with  respect  to  the  parameter  A  is 
denoted  by  the  superscript  '  ;  thus,  we  have  Xq'  =  dx0(A)/dA. 
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The  following  low-Froude-number  approximation  to  the  far-fteld  wave-amplitude 
function  K(t)  is  given  in  [3]: 

K(t)  ~  q2  (K+  +  K_)  as  F  -  0  ,  (4) 

where  the  port-and-starboard  contributions  K±  are  given  by 

K±  =  f  *s  exp  (-i6±/qF2)  a±  dA  .  (4a) 

ab 

In  these  equations,  F  is  the  Froude  number  defined  as 

F  =  U/(gL)i/2  ,  (5) 

and  we  have 

q  =  1/p  with  p  =  (l  +  t2)1/2  and  0  <  q  <  1 ;  (6a,b,c) 

furthermore,  0±  and  a±  are  the  phase  and  amplitude  functions  given  by 

0±  =  Xq  ±  ty0  and  (7) 

a±  =  u±a,±  +  F2q2(u±)2a2±  +  0(FA),  (8) 

where  u±  is  defined  as 

u±  =  l/[l-iq(x,±ty,)]  ,  (9) 

and  the  functions  a^  and  ^  are  now  defined. 

The  first-order  amplitude  function  a^  is  given  by 
*1*  =  yo'A±/(l  +  c2)  +  2q(x0'±ty0')(u±)2B±^0 

+  C±$0'  +  uD±^/(1  +  £2)  +  ip(y]V-yo'$i)*  (10) 

where  we  have 

£  =  (yo'xi-xo'yi)/11*  (H) 

with  u  defined  as 

u  =  KV)2  +  (yo')2l1/2>  (12> 

and  the  coefficients  A±,  B±,  C±  and  D±  are  given  by 
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(13a) 


A±  =  [(1 + py0  '/u)(l  -  py0  Vu) + 12]  +  Kpyo'/uJtyjCxo'ityo'yu  +  c] , 

B±  =  q(y2:?:tX2)  4  i(yix2-x1y2),  (13b) 

C±  =  [(Tt-p^'yo'/u2)  4  (py0  '/u)2e(x  ,  xq  ' + y  ,y0 '  )/u(l  +  c2)] 

+  i[y1(x0'±ty0')/u  +  £][px07u-(py07u)£(x1x0'+y1y0')/u(l-l-£2)] ,  (13c) 

Da:  =  [(x0'^yo')/u][(l  +  £2)(y1±iqt)u±+i(py07u)£y1]  -  (py07u)£(py07u-ic).  (13d) 

The  second-order  amplitude  function  >s  given  by 

a2±  =  — y  j '  +  2u±(<fom2±  +  $imi±  +  ^2mo±_iyo,y2±) 


4  6i(u±)2[^0(m1±y2±-l-m0±y3±)4  41m0±y2±l  "  ttfu^Wfo*)2  ,  (14) 

where  we  have 

Xn*  =  q^n^yn)  and  Mil*  =  Mn  +  iq(yn '  *Xn ' ) .  (15a, b) 

with 

Mo  =  Xiy0'-yix0'  =  at .  (16a) 

Mi  =  x^j'-yjx,'  +  2(x2y0'-y2x0'),  (16b) 

M2  =  xiy2-yix2'  +  2(x2y1'-y2x1')  4  3(x3y0'-y3x0') .  (16c) 

At  a  point  where  the  phase  of  the  trigonometric  function  exp(-i0±/qF2)  is 

stationary,  that  is  at  a  point  (xq  ,  y0 , 0)  where  we  have  Xo'±ty0'  =  0,  the  first-order 
amplitude  function  aj*  takes  the  form 

a,±  =  ±£[iy0'4(py1^itH0'Tu«fr1]/(l±i£)  4  ip(y^o'-yo'^i)  if  V^o'  =  °-  (17) 


Furthermore,  if  the  hull  intersects  the  plane  z  =  0  orthogonally  at  a  point  of 
stationary  phase  we  have 

a^  =  iF2py0'  d2<j>0/3x2  if  Xo'±ty0'  =  0  and  nz  =  0,  (18) 

and  the  amplitude  function  a±  then  is  of  order  F2. 

The  low-Froude-number  approximation  to  the  far-field  wave-amplitude  function 
K(t)  given  by  the  waterline  integral  (4)  is  well  suited  for  numerical  evaluation,  as  is 
shown  in  [3].  However,  expression  (4)  can  be  simplified  by  applying  the  method  of 
stationary  phase,  which  takes  advantage  of  the  rapid  oscillations  of  the  trigonometric 
function  exp(-i0±/ql'2)  in  the  low-Froude-number  limit  F  -*■  0.  General  presentations 


of  the  method  of  stationary  phase  may  be  found  in  Erdelyi  [4]  and  Copson  [5].  A 
brief  presentation  of  this  well-known  method  of  approximation  is  included  in  the 
following  section  for  easy  reference. 


THE  METHOD  OF  STATIONARY  PHASE 

Let  us  consider  the  integral 

1  =  f.  Xl  exp[iv0(t)]  a(t)  dt ,  (19) 

Jti 

where  v  is  a  large  positive  real  variable  and  the  phase-function  0(t)  is  real;  both  the 
phase-function  0(t)  and  the  amplitude-function  are  assumed  to  be  differentiable  in  the 
range  of  integration  [t(  ,  t2]  to  the  order  required  by  the  analysis.  We  seek  an 
asymptotic  expansion  of  the  integral  (19)  in  the  limit  v  — «.  Due  to  cancellations 
between  the  positive  and  negative  values  of  the  rapidly-oscillating  function  exp[iv0(t)], 
the  major  contributions  to  the  value  of  the  integral  I  stem  from  the  immediate 
vicinity  of  the  end  points  tj  and  t2  of  the  integration  range,  on  one  hand,  and  from 
the  vicinity  of  those  points  where  the  phase-function  0(t)  is  stationary,  that  is  where 
0'(t)  =  0,  on  the  other  hand. 

The  contribution  of  the  end  points  may  be  determined  by  supposing  that  the 
phase-function  0(t)  is  not  stationary  within  the  integration  range,  that  is  we  now 
assume  0'(t)  ¥=  0  for  tj  <  t  <  t2.  By  integrating  the  integral  (19)  by  parts,  we  may 
obtain 


I  =  _L  “<0-  eiv«,)  ‘2  _  ±  f<2  cM  ( aydt . 
iv  d'(t)  t,  iv  ■'t,  lev 

Another  integration  by  parts  yields  the  two-term  expansion 


l  oj  _L  eiv0(t) 

iv  0'(t) 


t2 . 1MI', i 


ivd 


|tj  V2  0' 
which  may  be  expressed  in  the  form 


*2 


I  ~ 


as  v 


00 


(20) 


If  the  functions  0(t)  and  a(t)  are  sufficiently  smooth  at  the  end  points,  higher-order 
terms  in  the  asymptotic  expansion  (20)  can  be  obtained  by  continuing  the  process  of 
integrating  by  parts. 

Let  us  now  assume  that  the  phase-function  0(t)  is  stationary  at  the  interior  point 
t  =  T,  where  tj  <  T  <  t2;  we  then  suppose  0'(T)  =  0,  but  0"(T)  #  0.  The  notation 

0  =  0(T),  0"  =  0"(T) ,  ...  A  =  a(T) ,  A'  =  a'(T),  ...  (21) 
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will  be  used  for  simplicity.  Let  the  function  <p(t)  be  defined  as  <p(t)  =  0(t)  -  0(T)  = 
6(t)  -  0.  By  performing  the  change  of  variable  t  =  t  -  T  we  may  then  express  the 
integral  I  in  the  form 

I  =  exp(iv0)  J,  (22) 

where  J  is  the  integral  defined  as 

J  =  f  2  exp[iv<p(T+r)]  a(T+T)  dT 
J  -Tj 

with  Tj  =  T  —  tj  and  t2  =  t2  -  T,  so  that  we  have  tj  >  0  and  t2  >  0.  By 
expanding  the  functions  <p(T+t)  and  a(T+r)  in  Taylor  series  about  the  origin  t 
we  may  obtain 

<p(T+r)  =  t20"/2  +  t3©"76  +  t4©(4>/24  +  ... 
and  a(T+T)  =  A  +  tA'  +  t2A"/2  +  ...  .  (24) 

The  change  of  variable  <p(T+t)  =  u20"/2  yields 

u2  =  t2  +  t30'"/30"  +  t4©^)/^©"  +  ...  . 

By  inverting  this  series  we  may  obtain 

t  =  u  -  c2u2/3  +  c3u3/3  +  ...  ,  (25) 

where  the  coefficients  c2  and  c3  are  given  by 

c2  =  ©"720"  and  c3  =  [5(0"')2/3©''-0(4>]/80"  .  (26a,b) 

By  performing  the  change  of  variable  <p(T+t)  =  u20'72  in  the  integral  (23)  we  may 
obtain 

J  =  f^2  exp[iv0''u2/2]  a(u)  du  ,  (27) 

7-u, 

where  the  function  a(u)  is  defined  as 

a(u)  =  3(T+t)  dr/du.  (28) 

By  using  the  Thyior  series  (24)  and  (25)  in  equation  (28)  we  may  obtain 

O’(u)  =  A  +  u(A'-2c2A/3)  +  u2( A' 72  -  c2A'  +  c3 A)  +  ...  .  (29) 

Equations  (27)  and  (29)  then  yield 

J  =  AJ0  +  (A'-2c2A/3)J!  +  (A'  72  -  c2A'  +  c3 A) J2  +  ...  ,  (30) 
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where  Jn  are  the  integrals  defined  as 

Jn  =  f  “2  exp[iv0"u2/2]  un  du  .  (31) 

J  -Uj 

The  change  of  variable  v0"u2/2  =  ct2,  where  t  is  defined  as 

t  =  sign  6"  ,  (32) 

yields 

Jn  =  (2/v|0"|)(n+O/2ln,  (33) 


where  the  integral  In  is  defined  as 

In  =  f *2  exp(ict2)  tn  dt 
J  -tj 

with  tj  =  (v|0"|/2)1/2Uj.  We  then  have  tj  -*  °°  as  v  -*■  oo.  The  integral  I0  is  given  by 

I0  ~  f  exp(ict2)  dt  =  n1/2  exp(i£rr/4) .  (34) 

The  integral  In  is  not  defined  in  the  limit  tj  -*  00  for  n  >  1.  However,  the  contribution 
of  the  immediate  vicinity  of  the  point  of  stationary  phase  t  *  T  to  the  integral  (19) 
can  be  evaluated  by  expressing  the  integral  In  in  the  form 

I n  =  f  exp(i£t2)  tn  o(t)  dt ,  (35) 

✓  —  00 


where  the  function  o(t)  is  equal  to  1  for  finite  values  of  t  but  vanishes  exponentially 
as  t  -*  By  performing  an  integration  by  parts  we  may  obtain 


21  j  =  -U  exp(Ut2)  o(t) 


+  i £  /  exp(iet2)  o'(t)  dt . 
—  00  •/  — -  00 


(36) 


This  yields 


I,  =  0 

since  we  have  o(t)  -*•  0  as  t  -* 
2I2  =  -U  exp(i£t2)  to(t) 
This  yields 


±  oo  and  o'(t)  =  0.  We  may  similarly  obtain 

oo  r  oo 

+  i £  /  exp(i£t2)  d[to(t)J . 

—  oo  •/  —  00 


2I2  = 


exp(i£t2)  o(t)  dt  =  iel0  , 


(37) 


(38) 


(39) 
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as  may  be  seen  from  Eq.  (34).  By  using  Eqs.  (34),  (37)  and  (39)  in  Eq.  (33)  we  may 
obtain 

J0  =  (2n/v|6"|),/2  exp(itn/4) ,  Jj  =  0,  J2  =  Uq/vQ"  .  (40a, b,c) 

Equations  (40a, b^),  (30)  and  (22)  then  yield 

I  ~  (2n/v|0"|),/2  exp(iv0+ieti/4)  {A+iA2/2v+0(l/v2)]  as,  v  -  »  ,  (41) 

where  A2  is  defined  as 

A2  =  (A"-2c2A,+2c3A)/0". 

Equations  (26a, b)  finally  yield 

A2  =  [(a70")'+a{5(0''')2/30"-e<4)}/4(e")2]t=x.  (42) 

The  contribution  of  an  interior  point  of  stationary  phase  t  =  T,  with  tj  <  T  <  t2, 
thus  is  given  by  Eq.  (41),  where  t  is  equal  to  the  sign  of  0"(T)  and  A2  is  defined  by 
Eq.  (42). 

An  analytical  approximation  to  the  integral  (19),  valid  in  the  limit  v  -*  °°,  is  then 
given  by  the  sum  of  the  asymptotic  expansion  (20),  corresponding  to  the  contribution 
of  the  end  points  tj  and  t2,  and  of  the  asymptotic  expansion  (41)  for  each  interior 
point  of  stationary  phase  t  =  T,  with  ti  <  T  <  t2.  Equations  (20)  and  (41)  show  that 
the  contributions  of  the  end  points  and  of  a  point  of  stationary  phase  are  0(1 /v)  and 
0(l/v1/2)  as  v  -*•  °°,  respectively.  The  major  contribution  to  the  integral  (19)  therefore 
stems  from  the  end  points  only  if  there  is  no  point  of  stationary  phase  within  the 
integration  range.  The  asymptotic  expansions  (20)  and  (41)  are  not  valid  in  the  special 
case  when  the  end  point  tj  or  t2  is  a  point  of  stationary  phase,  that  is  if  we  have 
0'(tj)  =  0  or  0'(t2)  =  0.  In  this  special  case,  the  two  asymptotic  expansions  should 
be  replaced  by  the  expansion  given  below. 

Let  us  suppose  that  the  phase-function  0(t)  is  stationary  at  the  end  point  tj  or  t2, 
so  that  we  have  0'(tj)  =  0,  but  0"(t,)  #  0.  The  contribution  of  this  end  point  of 
stationary  phase  may  easily  be  obtained  from  the  foregoing  analysis  for  an  interior 
point  of  stationary  phase.  More  precisely,  the  contribution  of  the  end  point  tj  is  given 
by  equations  (22),  (30),  (33)  and  (35)  where  the  range  of  integration  [  _  oo  ,  oo]  must 
simply  be  replaced  by  [0 ,  oo]  or  [  _  »  ,  0]  if  the  phase  is  stationary  at  the  lower  or 
upper  end  point  tj  or  t2,  respectively.  It  may  then  be  verified  that  the  values  of  the 
integrals  I0  and  I2  become  equal  to  half  the  values  given  by  Eqs.  (34)  and  (39),  and 
the  value  of  the  integral  Ij  is  equal  to  Tit/2  where  the  upper  sign  (-)  and  the  lower 
sign  (+)  correspond  to  the  cases  when  t2  and  tj  are  stationary  points,  respectively. 

The  contribution  of  the  end  point  tj  when  the  phase  is  stationary  there  is  then  given 
by 
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(43) 


I  ~  (n/2v|0  ‘|)1/2  cxp(iv©+icn/4) 


I 


[A  T  it  cxp(  - fcn/4) A, /(2v) 1  /2  +  iA2/2v  +  0(l/v3/2)]  as  v 


00 


< 

where  t  is  equal  to  the  sign  of  0"  at  the  end  point,  the  upper  sign  (-)  and  the  lower 
sign  (+)  correspond  to  the  cases  when  the  phase  is  stationary  at  the  upper  and  lower 
end  point,  respectively,  At  is  defined  as  j 


A!  =  2[(a'-a0'"/30")/(7i|0"|)1/2lt_t.  , 


(44) 


and  A2  is  given  by  Eq.  (42)  with  T  =  tj . 

The  asymptotic  expansion  (41)  for  the  contribution  of  an  interior  point  of 
stationary  phase  t  =  T  is  valid  if  the  first  derivative  0'(t)  of  the  phase  function  0(t) 
vanishes  at  the  point  t  =  T  but  the  second  derivative  0"(t)  is  nonzero  there,  as  was 
already  noted.  However,  the  asymptotic  expansion  (41)  is  clearly  not  valid  if  both  the 
first  and  the  second  derivatives  vanish  at  t  =  T.  An  asymptotic  expansion  for  the 
contribution  of  an  interior  point  of  stationary  phase  t  =  T,  with  tj  <  T  <  t2>  in  the 
case  when  0'(T)  =  0  and  0''(T)  =  0,  but  0"'(T)  #  0,  can  be  obtained  in  a  manner 
similar  to  that  used  for  obtaining  the  asymptotic  expansion  (41).  The  change  of 
variable  cp(T+  t)  =  u2©"/2  used  below  Eq.  (24)  and  prior  to  Eq.  (27)  is  replaced  by 
<p(T+t)  =  u30"'/6,  and  Eqs.  (25)  and  (26a, b)  become 


t  =  u  -  c2u2/6  +  c3u3/6  +  ...  with  (45) 

c2  =  eW/20'"  and  c3  =  [(0<4))2/40"'-0(5>/5)/20"' .  (46a, b) 


Equations  (30)  and  (31)  then  take  the  form 

J  =  AJ0  +  (A'-c2A/3)Ji  +  (A"-c2A'+c3A)J2/2  +  ... 
with  Jn  =  fUl  exp[iv0"'u3/6]  un  du  . 

•'-Uj 

The  change  of  variable  v0”'u3/6  =  t3  yields 
Jn  =  (6/v©"')(n+1>/3  In 

/OO 

exp(it3)  tno(t)  dt  . 

-  00 

It  may  be  verified  that  we  have 

I0  =  f(l/3)/31/2  ,  i!  =  i2n/3f(l/3)  and  I2  =  0 , 


(47) 

(48) 


(49) 

(50) 

(51a, b,c) 


where  T(l/3)  =  2.6789...  .  Equations  (22),  (47),  (49),  (51a, b^)  and  (46a, b)  finally  yield 
the  asymptotic  expansion 


9 


‘»  ut «««  “A tthf  A u  ftMRMAkJfUtAK  AM nxAJTAJCTDCTUfZUe  ^  TV*  -IVtt f\X  AJtTUTawJ 


I  ~  3-,/2r(l/3)(6/v0'"),/3  exp(ivQ)  [A  +  iA2/v1/3  +  0(l/v)]  as  (52) 

where  A2  is  defined  as 

A2  =  [ti24/3  /  3,/6{r(l/3)}2]  [(a a#4)/ 60 ' ' ' )/(0 " ' )  1/3]t  =T  .  (53) 

This  asymptotic  expansion  defines  the  contribution  of  an  interior  point  of  stationary 
phase  t  =  T,  with  tj  <  T  <  t2  ,  where  we  have  0'(T)  =  0  =  0"(T)  but  0"'(T)  #  0. 


BOW-AND-STERN  CONTRIBUTION 

Equations  (4a),  (19)  and  (20)  show  that  the  contribution  of  the  bow  and  the  stern 
to  the  low-Froude-number  asymptotic  approximation  to  the  integral  K±  is  given  by 

l*S 


K±  ~  (iF2q/0;)  [a±-iF2q(a±/e;)'  +  0(F4)]  exp(-ip0±/F2) 
By  using  Eqs.  (4),  (7)  and  (8)  we  may  then  obtain 
K  ~  iF2q3  [ABexp(  -  ipxB/F2)  - Asexp(  -  ipxs/F2)] , 


(54) 


where  x^s  are  the  abscissae  of  the  bow  and  the  stern  and  the  bow-and-stern 
amplitude  functions  A^s  are  defined  as 

Abs  =  2nxA,  +  iF2q(A2++A2-)  +  0(F4)  ,  (55) 

where  the  expression  on  the  right  side  is  evaluated  at  the  bow  or  at  the  stern.  In  this 
equation,  the  second-order  amplitude  functions  Af  are  given  by 

Af  e;  =  (u^/e;)'  +  iq(u±)V  .  (56) 


and  the  first-order  amplitude  function  Aj  is  defined  as 
-2nx  A!  =  u  +  ai+/0+  +  u_af/0'_ 


(57) 


for  reasons  that  will  become  clear  further  on. 

Equation  (2b)  yields  yj  =  0  since  we  have  y  =  0  at  the  bow  or  the  stern. 
Equation  (9)  thus  becomes 

u+  =  1/(1  —  iqxj)  =  u_  .  (58) 

Equation  (57)  then  yields 


-2nx(l-iqX|)0V®'- Aj  =  9'+  . 
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By  using  Eqs.  (7),  (10),  (58)  and  the  relation  yj  =  0  into  the  foregoing  equation,  we 
may  obtain 

-Inxd-iqx^eveLA,  =  yo'(A+0'_  +A_0;)/(l  +  £2) 

+  2qe;euB++B_)V(i-iqx,)2  +  (c+ei  +c_e;w 
+  u(D+0L  +D_eVH2/(l  +  £2)  -  2ipx<Sy^1 .  (59) 

By  using  the  relation  yj  =  0  into  Eqs.  (13a)  and  (11)  we  may  obtain 

A±  =  (1 + pyo/u)(l  -  py^u)  +  £2  +  icpyo/u  =  A  ,  (60) 

with  t  =  Xjyo/u.  (61) 

We  then  have 

a+ol  +a_ov  =  A(e;  +ei)  =  2x<^a  , 

where  Eq.  (7)  was  used.  By  using  Eqs.  (60)  and  (61)  we  may  finally  obtain 

A+flL+A_e;  =  2x^[l+(x,2-p2+ipx,)(y6/u)2] .  (62) 

By  using  the  relations  y!  =  0  and  y2  =  0,  which  follow  from  Eq.  (2b)  and  the 
condition  y  ==  0  at  the  bow  or  the  stern,  in  Eq.  (13b)  we  may  obtain  B±  =  Tqtx2.  We 
then  have 

B+  +  B_  =0.  (63) 

By  using  the  relation  yj  -  0  into  Eq.  (13c)  we  may  obtain 
C±  =  Tt  -  (py,)/u  -  ie)[px^/u  -  (py^/u)£x1xyu(l  +  £2)]  , 
which  becomes 

C±  =  Tt  -  (l-iqx,)(p2x^y(J/u2)/(l  +  £2) 
by  virtue  of  Eq.  (61).  By  using  Eq.  (7)  we  may  finally  obtain 

C  +  01  +C_e;  =  2yg[t2 -(1  -iqx,)(pxQ/u)2/(l  +  £2)] .  (64) 

Equations  (13d),  (58),  (7)  and  the  relation  yj  =  0  yield 
D±  =  ±iqt(l  +  £2)©i/u(l  - iqxj)  -  (py^/uMpy^/u-if) . 


We  then  have 


D  +  61+D_0V  =  -2p2x1x^(y^/u)3(l-iqx1),  (65) 

where  Eqs.  (7)  and  (61)  were  used. 

By  substituting  Eqs.  (62),  (63),  (64)  and  (65)  into  Eq.  (59)  we  may  obtain 

(l-iqxj)(l  +  c2),/2[(x^/u)2  -(ty<J/u)2]A,  =  -  [1 +(x,2  -  p2 +ipx1)(y(J/u)2]x(J/u 
-  [t2 + (tx i yo/u)2  -  (1  -  iqx  j  )(px^/ u)2] u  +  ip(l — ipx,(y^/u)2J41x^/u ,  (66) 

where  Eqs.  (7)  and  (61)  and  the  relation  nx  =  yo/u(l  +  £2)1/2,  given  by  Eq.  (33a)  in 
(3],  were  used.  Equations  (31)  and  the  equation  below  it  in  [3]  yield 

xjtx+t  -(txny -tynx  +x,tynz)$d  ,  (67) 

where  the  relation  yj  =  0  was  used.  Furthermore,  we  have 

£  =  -nz/(l-nz2)1/2  and  l  +  £2  =  1/(1  -nz2)  ,  (68a, b) 

as  may  be  obtained  from  Eq.  (33c)  in  [3].  Equations  (39),  (38a, b)  and  (33a, b)  in  [3] 
also  yield 

-46/u  =  *t  ,  (69a) 

-Xq/u  =  tx  =  ny/(l-nz2)1/2  ,  (69b) 

-yo/u  =  ty  =  -nx/(l-nz2)1/2  ,  (69c) 

where  Eq.  (68b)  was  used.  Equations  (61)  and  (69c)  give 

Xj  =  -nz/nx  .  (70) 

Equation  (67)  becomes 

♦j  =  -txnzVnx-V0-nz2)1/2  (71) 

upon  using  Eqs.  (70)  and  (69b, c).  By  substituting  Eqs.  (69a, b,c),  (70)  and  (71)  into  Eq. 
(66)  we  may  obtain 

(nx  +  iqnz)(tx2-t2ty2)A,  =  ny(nx2+ty2(nz2-p2nx2-ipnxnz)]/nx 

-  ty(t2  -  p2tx%  +  iptx(nx + ipty2nz)f d  . 

By  using  Eqs.  (6b),  (69b, c)  and  the  relations  tx2  +  ty2  =  1  and  nx2  +  ny2  +  nz2  =  1  in 
the  foregoing  expression  we  may  obtain 
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(nx+iqnz)A!  =  nxny  [1  -  ipn2  (nx  +  iqnz)/(l  -  nz2  -  p2nx2)] 

+  tyjt  -  ipnyty*d[l  +  p2nx(nx  +  iqrg/(l  -  nz2  -  p2^2)]  . 

Upon  dividing  by  nx+iqnz  we  may  finally  obtain 

A,  =  nx(ny  -  vt  +  ipnyVd)/(nx + iqnz) 

+  ipnxny(nxipd  -  q2nz)/[q2(l  -  nz2)  -  nx2]  ,  (72) 

where  we  have 

(V,  ,  Vd>  =  rtt  •  td)/U-nz2)1/2  •  (73a, b) 

The  contribution  of  the  bow  and  the  stern  to  the  far-field  wave-amplitude 
function  K(t)  is  then  given  by  Eqs.  (54)  and  (55),  where  the  first-  and  second-order 
amplitude  functions  A,  and  A2±  are  defined  by  Eqs.  (72)  and  (56),  respectively. 
Expression  (56)  for  the  second-order  amplitude  function  is  complex.  However,  Eq. 

(72)  defines  the  first-order  amplitude  function  Aj  explicitly  in  terms  of  t,  p  = 

(l+t2)1/2  and  q  =  1/p  ,  the  geometrical  characteristics  of  the  hull  at  the  bow  and  the 
stern  —  specifically,  the  unit  outward  hull-normal  vector 7f(nx,ny,nz)  —  and  the 
components  +t  and  <{»d  of  the  disturbance  velocity  vector  in  the  directions  of  the  unit 
hull-tangent  vectors-T and  nxT 

Equation  (55)  shows  that  the  first-order  approximation  to  the  bow-or-stern 
amplitude  function  A^s  —  given  by  2nxAj  —  vanishes  if  nx  =  0,  that  is  if  the 
waterline  has  a  cusp  at  the  bow  or  the  stern.  Equation  (72)  shows  that  the  first-order 
approximation  to  the  amplitude  function  A^s  vanishes  also  if  the  bow  or  the  stem  is 
round,  since  we  then  have  ny  =  0  and  $t  =  -  4>y  =  0  by  symmetry.  It  may  thus  be 
seen  that  the  contribution  of  the  bow  (or  the  stern)  to  the  far-field  wave-amplitude 
function  K(t)  is  of  order  F2  if  nx  =  0  or  ny  =  0  at  the  bow  (or  the  stem);  that  is,  we 
have 

Kas  =  0(F2)  if  nxB>s  =  0  or  nyB-s  =  0  .  (74) 

In  the  particular  case  when  the  hull  surface  is  vertical  at  the  bow  or  the  stern,  we 
have  nz  =  0  and  Eq.  (72)  becomes 

Ai  =  lx  ~  ~  iqtxVfa2-^2)  if  nz  =  0,  (75) 

where  the  relations  <{>d  =  -  and  ny  =  tx  were  used.  This  yields 

A,  =  tx  -  *t  +  0(F2)  if  nz  =  0 ,  (76) 

as  is  indicated  by  the  free-surface  boundary  condition  <j>z  =  -F2^  =  0(F2). 

In  the  limit  t  -*•<»,  we  have  q  -*  0  and  Eq.  (72)  yields 
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'  -*  ny  -  (<t>t  -  <|>dnzny/nx)/ (1  -  nz2) 1/2  as  t 


00 


(77) 


lu  the  particular  case  t  =  0,  we  have  p  =  1  =  q  and  Eq.  (72)  becomes 

A,  =  nx(ny-ipt+inyipd)/(nx+iqnz)  +  inx(nxi/;d-nz)/ny  for  t  =  0  .  (78) 

Expression  (78)  is  not  valid  if  ny  =  0,  that  is  if  the  bow  or  the  stern  is  round. 
More  generally,  expression  (72)  is  not  valid  if  q2  =  nx2/(l-nz2)  =  ty2  ,  that  is  if 
tx2  -  t2ty2  =  0.  This  special  case,  in  which  the  bow  or  the  stern  is  a  point  of 
stationary  phase,  is  considered  further  on. 

STATIONARY-PHASE  CONTRIBUTION 

Equations  (4),  (4a),  (7),  (8),  (19),  (41)  and  (42)  show  that  the  contribution  of  the 
point(s)  of  stationary  phase  may  be  expressed  in  the  form 

K  ~  i  F  q3  E  K±  ;  (79) 

in  this  expression,  K±  corresponds  to  the  contribution  of  a  point  of  the  mean 
waterline  where  the  phase  of  the  trigonometric  function  exp[-ip(xo±ty0)/F2]  is 
stationary,  that  is  where  we  have 

t  =  ^/dyo  =  **<^6  =  Ttx/ty  ,  (80) 

as  may  be  obtained  by  using  Eqs.  (69b, c),  and  the  notation  2-  implies  summation 
over  all  the  points  of  stationary  phase. 

The  stationary-phase  contribution  K±  in  Eq.  (79)  is  given  by 

K±  =  ±(2nr),/2  A±  exp[-i(p0±/F2T£7i/4)]  ,  (81) 

where  r  and  z  are  defined  as 


*- 


r  =  u2/q|xo'±tyo  |  , 
z  =  t  sign  (xo'±tyo')  ; 

furthermore,  the  amplitude-function  A±  is  given  by 


(82) 

(83) 


(86) 


±2uA2 *  =  Ku^j^'/e;']'  +  U±a ,^[5(04' ' >2/304'-  /4<e±  )2 

+  2iq(u±)2a2±  . 

In  Eqs.  (82),  (85)  and  (86),  we  have  u  =  [(xq)2  -t-  (yo)2] 1/2.  as  is  given  by  Eq.  (12); 
furthermore,  we  have  6±  =  xq±  ty0,  as  is  indicated  by  Eq.  (7),  in  Eqs.  (81)  and  (86). 
Equations  (6a, b),  (12),  (69c)  and  the  stationary-phase  relation  (80)  yield 

q  =  ±y<>/u  =  Tty  =  1/p,  (87) 

where  the  condition  <  0,  which  follows  from  Eq.  (la),  was  used.  By  using  Eqs.  (7), 
(80)  and  (87)  we  may  then  obtain 

pfl*  =  ±(y0tx-xoty)/ty2  and  (88) 

r  =  [(xo)2  +  (yo)2]3/2/|xo  yo  -yo'xol  .  (89) 

which  shows  that  r  is  equal  to  the  radius  of  curvature  of  the  mean  waterline  at  the 
stationary-phase  point  (xq  ,  y0). 

Let  (a  ,  ft )  be  the  (x  ,  y)  coordinates  of  the  center  of  curvature  of  the  mean 
waterline  at  the  point  of  stationary  phase  ,  y0).  We  have 

*  =  Xo+y6[(x6)2+(y6)2]/(x6'y6-yo'x6)  • 

Equations  (83)  and  (80)  yield 

t  t  sign  (x^y<i-y6'x6)/y6  • 

We  then  have 

£  =  =F  sign  (a-xo)  ,  (90) 

which  shows  that  £  takes  the  values  -1  or  +1  if  the  center  of  curvature 
corresponding  to  the  stationary-phase  point  (xo  ,  y0)  is  upstream  or  downstream  from 
(xo  ,  yo),  respectively.  Equation  (80)  shows  that  we  have  ?  ty  ^  0  at  a  point  of 
stationary  phase,  since  we  have  t  >  0  and  tx  >  0.  More  precisely,  Eq.  (7)  shows  that 
the  points  where  the  phases  6+  and  6_  are  stationary  are  located  in  the  fore  and  aft 
portions  of  the  mean  waterline,  where  we  have  ty  <  0  and  ty  >  0,  respectively.  If  the 
mean  waterline  is  convex,  we  have  a  <  Xq  and  a  >  in  its  fore  and  aft  portions, 
respectively.  It  may  then  be  seen  from  Eq.  (90)  that  t  is  equal  to  +1  or  -1  if  the 
mean  waterline  is  convex  or  concave,  respectively,  at  the  point  of  stationary  phase. 

The  relation  -x^/u  =  [l-(y6/u)2]1/2.  which  follows  from  Eq.  (12),  and  Eqs.  (87) 
and  (6a,b)  yield 

qt  =  -x<j/u  =  tx  ,  (91) 
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where  Eq.  (69b)  was  used.  By  using  Eqs.  (87)  and  (91)  in  Eq.  (11)  we  may  obtain 


c  =  ±q(x!±tyj) .  (92) 

Equation  (9)  then  becomes 

u±  =  1/(1  Tk) .  (93) 

Upon  using  Eqs.  (17),  (93),  (69a, b/:)  and  (87),  Eq.  (85)  becomes 

(l  +  c^Aj*  =  ±qc  ±  (te— py,)4t  -  .  (94) 

Equation  (31)  and  the  equation  below  it  in  [3]  yield 

+1  =  (^x+yityWt  “  D+d  •  (95) 

where  D  is  defined  as 

D  =  txny-tynx+(x,ty-y1tx)n2  .  (96) 

Equations  (69b, c),  the  relation  tx2  +  ty2  =  1  and  Eqs.  (11)  and  (68a, b)  show  that  we 
have 

D  -  (l  +  £2)l/2  .  (97) 

Upon  using  Eqs.  (97),  (87)  and  (91)  in  Eq.  (95)  we  may  obtain 

+1  =  qftx^iWt  -  (l  +  c2),/2$d  .  (98) 

It  follows  from  Eqs.  (2a, b)  that  the  vector  d~x/dz  =  (xj  ,  yj  ,  1)  is  tangent  to  the 
hull  surface.  We  then  have  rT-  dx/dz  =  0,  that  is 

x,nx  +  yjny  +  nz  =  0  .  (99) 

This  equation  yields 

X,  =  Ttyj  -  nz/nx  ,  (100) 


where  Eqs.  (69b, c)  and  (80)  were  used.  By  using  Eqs.  (100)  and  (6a, b)  in  Eq.  (98)  we 
may  obtain 


-<h  =  (qtnz/nx±py,)<frt  +  (l  +  £2),/2$d  . 


(101) 


Equations  (94)  and  (101)  then  yield 

(l  +  c2)Ai*  =  ±qt  +  t(qnz/nx±c)4t  +  (l  +  t2),/2fd  .  (102) 

Equations  (68a),  (69c)  and  (87)  show  that  we  have 
qnz/nx  ±  t  =  0  . 

Equation  (102)  finally  becomes 

A,±  =  -nznx  +  (l-nz2)1/24d  ,  (103) 

where  Eqs.  (68a, b),  (69c)  and  (87)  were  used. 

Upon  using  Eq.  (88)  we  may  express  Eq.  (81)  in  the  form 

K±  =  ±(2nr),/2  A±  exp[Ti{(y0tx  -Xoty)/ty2F2  -  rn/4}]  ;  (104) 

in  this  expression  for  the  contribution  of  a  point  of  stationary  phase  (xq  ,  yo)»  which 
is  defined  by  Eq.  (80),  r  represents  the  radius  of  curvature  of  the  mean  waterline  at 
the  point  (xo  ,  yo),  c  is  equal  to  +1  or  -1  if  the  mean  waterline  is  convex  or  concave, 
respectively,  at  (xq  ,  y0)  and  A±  is  the  amplitude  function  defined  by  Eqs.  (84),  (86) 
and  (103). 

Expression  (86)  for  the  second-order  amplitude  function  A 2±  is  a  complex  one. 
However,  Eq.  (103)  defines  the  first-order  amplitude  function  Aj±  explicitly  in  terms 
of  the  geometrical  characteristics  of  the  hull  —  specifically,  the  components  nx  and  nz 
of  the  unit  outward  normal  vector  tT  to  the  hull  —  and  the  velocity  component  4d  in 
the  direction  of  the  downward  tangential  unit  vector  In  xT  at  the  point  of  stationary 
phase. 

Equations  (84)  and  (103)  show  that  in  the  particular  case  when  the  hull  surface  is 
vertical  at  the  point  of  stationary  phase,  we  have 

A±  =  -*z  +  0(F2)  if  nz  =  0.  (105) 

The  free-surface  condition  =  -F2^  =  0(F2)  then  yields 


A*  =  0(F2)  if  n2  =  0 . 


(106) 


The  summation  in  expression  (79)  for  the  contribution  of  the  point(s)  of 
stationary  phase  is  extended  to  any  point  of  the  mean  waterline  where  the  phases  of 
the  trigonometric  functions  exp[-ip(xo±tyo)/F2]  are  stationary,  that  is  where  Eq.  (80) 
holds.  The  number  of  stationary  points,  and  their  positions  along  the  waterline, 
depend  on  the  value  of  t  and  on  the  shape  of  the  waterline.  For  instance,  for  the 
simple  case  of  a  hull  with  waterline  consisting  of  a  sharp-ended  parabolic  bow  region 
1/4  <  x  <  1/2  defined  by  the  equation  y  =  4bx(l-2x),  where  b  denotes  the  ship’s 
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beam/length  ratio,  a  straight  parallel  midbody  region  -1/4  <  x  <  1/4,  and  a  round- 
ended  elliptic  stern  region  -1/2  <  x  <  -1/4  defined  by  the  equation 
y  =  b[-2x(l+2x)]1/2,  there  is  one  point  of  stationary  phase  in  the  stern  region  given 
by  x  =  -[l  +  l/(l+4b2t2)1/2]/4,  so  that  we  have  -1/2  <  x  <  -1/4  for  0  <  t  <  00 
with  x  -*•  — 1/2  as  t  — ►  0  and  x  -*  -1/4  as  t  and  one  point  of  stationary  phase 
in  the  bow  region  given  by  x  =  (1  +  l/4bt)/4  for  l/4b  <  t  <  °°,  so  that  we  have  1/2 
>  x  >  1/4  for  l/4b  <  t  <  00  with  x  -*  1/2  as  t  -*  l/4b  and  x  -*■  1/4  as  t  -*•  ».  We 
thus  have  one  point  of  stationary  phase  in  the  stern  region  for  0  <  t  <  l/4b  and  two 
points  of  stationary  phase,  one  in  the  stern  region  and  one  in  the  bow  region,  for 
l/4b  <  t  <  00 .  The  two  points  of  stationary  phase  approach  the  shoulders  x  =  ±1/4, 
where  dy/dx  =  0,  in  the  limit  t  -*  «>. 

Expression  (104)  is  not  valid  at  a  point  of  stationary  phase  of  order  two,  for 
which  both  the  first  and  the  second  derivatives  of  the  phase  vanish.  In  this  case,  we 
thus  have 

xq  ±  tyB  =  0  and  Xq'  ±  tyg  =  0  . 

These  two  conditions  can  be  satisfied  simultaneously  at  a  point  (xo  ,  y0)  where  we 
have 


x6'y6-y6,x6  =  0 . 

that  is  at  a  point  of  inflexion  of  the  mean  waterline,  at  which  the  radius  of  curvature 
is  infinite  as  may  be  seen  from  Eqs.  (89)  or  (82).  This  particular  case  is  considered 
further  on. 


SUMMARY  OF  RESULTS 

Equations  (79)  and  (54)  yield  the  following  low-Froude-number  asymptotic 
expansion  for  the  far-field  wave-amplitude  function  K(t): 

K(t)  ~  iFq3  [  Z  K±  +  F(Kb-Ks)  ]  as  F  -  0  ,  (107) 

where  F  is  the  Froude^number  defined  by  Eq.  (5),  that  is: 

F  =  U/(gL),/2  ,  (108) 

q  is  given  by  Eqs.  (6a,b,c),  i.e.: 

q  =  1/p  with  p  =  (l+t2),/2  and  0<  q  <  1  ,  (109a, b,c) 

and  K±)  KB  and  Ks  represent  the  contributions  of  the  point(s)  of  stationary  phase  on 
the  mean  waterline,  and  of  the  bow  and  the  stern,  respectively. 

The  notation  Z  implies  summation  over  all  the  points  of  the  mean  waterline 
where  the  phases  o^the  trigonometric  functions  exp[-ip(x±ty)/F2]  are  stationary,  and 
the  terms  K±  correspond  to  the  contributions  of  the  stationary  points  of  the  phases 
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x±ty,  respectively.  These  points  of  stationary  phase  are  defined  by  the  equivalent 
relations  (80),  that  is: 

dy/dx  =  ty/tx  =  Tl/t  ,  (110) 

where  tx  and  ty  are  the  components  of  the  unit  vectorT tangent  to  the  mean  waterline 
and  pointing  towards  the  bow.  The  number  of  points  of  stationary  phase,  and  their 
positions  on  the  mean  waterline,  depend  on  the  value  of  t  and  on  the  shape  of  the 
waterline.  We  usually  have  one  or  two  points  of  stationary  phase  for  typical  hull 
forms,  as  was  shown  in  the  previous  section. 

The  stationary-phase  contribution  K±  is  given  by  Eq.  (104),  that  is: 

K±  =  ±(2nr)1/2  A±  exp[Ti{(ytx-xty)/ty2F2-£Ti/4}]  ;  (111) 

in  this  expression,  r  represents  the  radius  of  curvature  of  the  mean  waterline  at  the 
point  of  stationary  phase  (x  ,  y),  t  is  equal  to  +1  or  -1  if  the  mean  waterline  is 
convex  or  concave,  respectively,  at  (x  ,  y)  and  A±  is  the  amplitude  function  defined  by 
Eq.  (84),  that  is: 

A±  =  A,  -  F2qA2±  +  0(F4)  .  (112) 

The  first-order  amplitude  function  Aj  is  given  by  Eq.  (103),  i.e.: 

A,  =  -nznx  +  (l-nz2)l/2+d  ,  (113) 

where  "n(nx  ,  ny  ,  nz)  is  the  unit  outward  normal  vector  to  the  hull  and  <{>d  is  the 
velocity  component  in  the  direction  of  the  downward  tangential  vector  nxT to  the 
hull.  The  second-order  amplitude  function  \2±  *n  Eq.  (112)  is  given  by  a  complex 
expression,  namely  Eq.  (86),  where  0±  ,  u±  ,  u  ,  aj*  and  a2±  are  defined  by  Eqs.  (7), 
(9),  (12),  (17)  and  (14).  Equations  (112)  and  (113)  show  that  in  the  particular  case 
when  the  hull  surface  is  vertical  at  the  point  of  stationary  phase,  we  have 

A±  =  -<j>z  +  0(F2)  =  0(F2)  if  nz  =  0  .  (114) 

The  bow  and  stem  contributions  KB  s  in  the  asymptotic  approximation  (107)  are 
defined  by  Eq.  (54),  that  is: 

kB,S  =  ab,S  exp(-ip  xb  s/F2)  ,  (115) 

where  xBS  are  the  abscissae  of  the  bow  and  the  stern,  respectively,  and  the  amplitude 
functions  ABS  are  defined  by  Eq.  (55),  i.e.: 

ABS  =  2nxAj  +  iF2q(A2+ +A2-)  +  0(F4)  .  (116) 
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The  first-order  amplitude  function  A]  is  given  by  Eq.  (72),  that  is 

Al  *  nx(ny  -  +  ipnyVd)/(nx + iqnz)  +  ipnxny(nx^d  - q2nz)/[q2(l  -  nz2)  - nx2]  ,  (117) 

where  (y>,  ,  yd)  =  ($t  ,  $d)/(l-nz2)1/2  (118a,b) 

and  $t  is  the  velocity  component  in  the  direction  of  the  tangent  vectorT to  the  mean 
waterline.  The  second-order  amplitude  functions  A2±  in  Eq.  (116)  are  given  by  a 
complex  expression,  namely  Eq.  (56),  where  9±  ,  u±  ,  a^  and  a2*  are  defined  by  Eqs. 
(7),  (9),  (10)  and  (14).  Equations  (116)  and  (117)  show  that  we  have  nxAj  =  0  and 
therefore 

ab,s  =  0(F2)  nx8,5  =  0  or  nyB>s  =  0 ,  (119) 

that  is,  if  the  bow  or  the  stem  is  either  cusped  or  round. 

The  asymptotic  approximation  (107)  and  Eqs.  (111)-(114)  and  (115)-(119)  defining 
the  contributions  of  the  point(s)  of  stationary  phase  on  the  waterline  and  of  the  bow 
and  the  stern,  respectively,  show  that  the  low-Froude-number  behavior  of  the  far-field 
wave-amplitude  function  is  strongly  influenced  by  the  shape  of  the  hull  in  the  vicinity 
of  the  mean  free  surface.  More  precisely,  for  a  value  of  t  for  which  there  is  one  (or 
more)  point  of  stationary  phase  on  the  mean  waterline  where  the  hull  has  flare,  the 
contribution  of  this  stationary-phase  point  dominates  the  contributions  of  the  bow 
and  the  stern  and  we  have  K(t)  =  O(F).  On  the  other  hand,  for  a  value  of  t  for 
which  either  there  corresponds  no  point  of  stationary  phase  or  the  hull  has  no  flare 
(i.e.  is  vertical)  at  the  point(s)  of  stationary  phase,  the  contributions  of  the  bow  and 
the  stern  are  dominant,  and  we  have  K(t)  =  0(F2),  except  if  both  the  bow  and  the 
stern  are  either  cusped  or  round.  In  the  latter  case,  the  contribution  of  the  point(s)  of 
stationary  phase,  where  the  hull  is  assumed  to  be  vertical,  is  dominant  and  we  have 
K(t)  =  0(F3). 

The  wave  resistance  of  a  ship  is  defined  in  terms  of  the  far-field  wave-amplitude 
function  by  means  of  the  well-known  Havelock  integral 

nR/CU2L2  «/"  ilK(t)ll2p  dt  .  (120) 

The  low-Froude-number  asymptotic  approximation  (107)  then  shows  that  the 
nondimensional  wave-resistance  coefficient  is  0(F2)  for  a  ship  form  with  a  region  of 
flare,  0(F4)  for  a  ship  form  that  is  vertical  (wall  sided)  everywhere  along  its  mean 
waterline  and  has  either  a  bow  or  a  stern  (or  both)  that  is  neither  cusped  nor  round, 
and  0(F6)  for  an  everywhere  wall-sided  ship  form  with  both  bow  and  stern  that  are 
either  cusped  or  round. 

In  the  case  of  a  ship  form  that  has  flare  over  a  portion  of  its  waterline  and  is 
wall  sided  elsewhere,  the  far-field  wave-amplitude  function  K(t)  is  0(F)  for  the  range 
of  values  of  t  for  which  there  corresponds  a  point  of  stationary  phase  within  the 
region  of  flare  but  K(t)  is  0(F2)  for  other  values  of  t,  for  which  the  corresponding 
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points  of  stationary  phase  fall  outside  the  region  of  flare.  If  the  portion  of  the 
waterline  where  the  ship  has  flare  is  of  very  small  extent,  the  corresponding  range  of 
values  of  t  for  which  the  function  K(t)  is  O(F)  is  also  small,  so  that  the  function  K(t) 
may  exhibit  a  sharp  peak  in  the  limit  F  -*■  0.  The  amplitudes  of  the  transverse  and 
divergent  waves  in  the  far-field  Kelvin  wake  of  a  ship  are  directly  related  to  the  far- 
field  wave-amplitude  function  K(t)  and  therefore  may  also  exhibit  a  sharp  peak  at  an 
angle  from  the  track  of  the  ship  smaller  than  the  Kelvin  cusp  angle  of  19°  1/2,  as  is 
shown  in  Noblesse  [2]. 


SPECIAL  CASES 

As  was  already  noted,  the  asymptotic  expansion  (107)  is  not  uniformly  valid  for 
all  values  of  t.  In  particular,  this  expansion  is  not  valid  in  the  vicinity  of  the  values  of 
t  for  which  the  bow  or  the  stern  is  a  point  of  stationary  phase,  that  is  for  t  =  tB  or 
ts  ,  with 

t[j  =  —  (dx/dy)jjOW  =  —  (tx/ty)bow  *  (121a) 

tg  =  (dx/dy)stern  =  (tx^y)stern  •  (121b) 

Equations  (41)  and  (43)  show  that  the  stationary-phase  and  bow  contributions 
K+  +  FKB  in  Eq.  (107)  become  K  +  /2  +  O(F)  for  the  special  value  t  =  tB.  The 
stationary-phase  and  stern  contributions  K  _  -  FKS  in  Eq.  (107)  likewise  become 
K_/2  +  O(F)  for  t  =  ts.  We  thus  have 

K+  +  FKb-  K  +  /2  +  O(F)  if  t  =  tB  ,  (122a) 

K_  -  FKS  -  K_/2  +  0(F)  if  t  =  ts  ,  (122b) 

where  K±  is  given  by  Eqs.  (1 11)-(1 13). 

Equations  (122a, b)  and  (1U)-(113)  yield 

K±  ±  FKbs  -*•  ±(nr/2)l/2  A  exp[±i(x/tyF2  +  7i/4)]  +  O(F)  if  t  =  tBS  ,  (123) 

where  y  and  t  in  Eq.  (Ill)  were  taken  equal  to  0  and  1,  respectively  (since  the 
waterline  must  be  convex  at  the  bow  and  at  the  stern),  and  the  amplitude  function  A 
is  given  by 

A  -  -nznx  +  (1  —  nz2)1/2  <t»d  .  (124) 

As  was  noted  previously,  the  asymptotic  expansion  (107)  is  also  not  valid  in  the 
vicinity  of  a  point  of  stationary  phase  for  which  the  mean  waterline  has  an  inflexion 
point;  this  corresponds  to  the  case  when  both  the  first  derivative  and  the  second 
derivative  of  the  phase  vanish,  that  is  we  have  both  Xq  ±  tyg  =  0  and  Xq'  ±  tyo'  =  0, 
which  yields  Xg'yQ  -  yo'xg  =  0,  i.e.  r  =  0  in  Eq.  (111). 
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An  expression  for  the  stationary-phase  contribution  K±  in  Eq.  (107)  valid  in  this 
special  case  can  be  obtained  from  Eq.  (52).  Specifically,  expression  (111)  for  the 
stationary-phase  contribution  becomes 

K±  =  ±F"1/3  3-1/2r(l/3)(6e)1/3  A±  exp[Ti(ytx -xty)/ty2F2]  ,  (125) 

where  q  is  defined  as 

q  =  [(x')2  +  (y')2]5/2/[y'(y'x'"-x'y'")J  ,  (126) 

and  the  amplitude  function  A±  is  given  by 

A±  =  Aj  +  0(F2/3)  ,  (127) 

with  Aj  given  by  Eq.  (113). 

CONTRIBUTION  OF  THE  NONLINEAR  TERMS  IN  THE 
FREE-SURFACE  BOUNDARY  CONDITION 

The  far-field  wave-amplitude  function  K(t)  considered  in  the  foregoing 
corresponds  to  the  usual  Neumann-Kelvin  theory,  in  which  the  free-surface  boundary 
condition  is  linearized.  The  far-field  wave-amplitude  function,  K'(t)  say,  associated 
with  the  generalized  Neumann-Kelvin  theory  in  which  the  nonlinear  terms  in  the  free- 
surface  boundary  condition  are  taken  into  account  is  defined  in  Noblesse  [1]  by  an 
expression  of  the  form 

K'(t)  «  K(t)  +  k(t)  ,  (128) 

where  k(t)  corresponds  to  the  nonlinear  terms  in  the  free-surface  boundary  condition. 
More  precisely,  Eqs.  (36)  and  (2)  in  [1]  show  that,  for  a  ship  with  port  and  starboard 
symmetry,  the  function  k(t)  is  given  by 

k(t)  =  k+(t)  +  k_(t)  with  (129) 

k±(t)  =  //f  exp(-i(x±ty)/qF2]  xW  dxdy  ,  (130) 

where  f  represents  the  positive  half  of  the  mean  free  surface  and  x(4)  represents  the 
nonlinear  terms  in  the  free-surface  boundary  condition  given  by 

x(4)  =  ia+/ax-(v*)2/2]  ao*/3z+F2a2*/ax2)/az 

-  a(v$)2/ax  +  vf V(v+)2/2  +  ckf2*3).  (i3i) 

Let  b  represent  the  beam/length  ratio  of  the  mean  wetted  hull  and  let  us  assume 
for  simplicity  that  a  line  parallel  to  the  x  axis  in  the  mean  free-surface  plane  z  =  0, 
defined  by  the  equation  y  =  rj  with  0  <  rj  <  b/2,  intersects  the  mean  waterline  at 
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only  two  points,  defined  by  x  =  4S  in  the  stern  region  and  x  =  |B  in  the  bow  region, 
say  (the  case  when  a  line  y  =  rj  intersects  the  mean  waterline  at  more  than  two  points 
yields  a  result  identical  to  that  obtained  below  and  may  be  treated  similarly). 

Equation  (130)  then  yields 


k±(t)  =yjj  dy  exp(?ity/qF2)  [  exp(-ix/qF2)  *(<}>)  dx 

+ /;  exp(  -  ix/qF2)  y(<j>)dx  j 

y*  +  oo  r  oo 

dy  exp(+ity/qF2)  J  ^  exp(  -  ix/qF2)  y(+)  dx  ■ 


The  inner  integrals  (with  respect  to  x)  in  the  foregoing  expression  are  Fourier 
integrals  for  which  asymptotic  expansions,  valid  in  the  limit  qF2  -*  0,  can  easily  be 
obtained,  e.g.,  from  Eqs.  (19)  and  (20)  or  from  Eqs.  (2)  and  (3)  p.  47  in  Erdelyi  [4], 
provided  that  the  free-surface  nonlinear  term  x(4)  and  its  x  derivatives  vanish  as  |x|  -* 
oo,  which  is  presumed.  We  may  then  obtain 

k±(t)  ~  iF2q[^b/2  exp[-i(|s±ty)/qF2]  [x  -  iF2q*x  +  0(F4)]^  dy 

-.T  exp[-i(|B±ty)/qF2]  [x-iF2qXx  +  °(f4)14b  dy  J  , 
which  may  be  expressed  in  the  form  of  the  line  integral 


/•  stern 

bow  exp[-i(x±ty)/qF2]  [x  - iF2qXx  +  0(F4)]  tydl  , 


(132) 


where  the  relation  dy  =  -tydl  was  used. 

By  using  the  parametric  representation  (la,b),  which  yields  tydl  =  -yBdA,  we  may 
express  the  waterline  integral  (132)  in  the  form 

k±(t)  ~  -iF2q  exp[-i(x0±ty0)/qF2]  [X - iF2qxx  +  0(F4)1  YodA  .  (133) 

J  Ab 

Equations  (128)  and  (4)-(8)  then  show  that  the  far-field  wave-amplitude  function  K'(t) 
in  the  generalized  Neumann-Kelvin  theory  may  be  expressed  in  the  form 


K'(t)  ~  q2(K'+  +  KL)  as  F—0  with  (134) 

K;(t)  ~  f^s  exp[  -  i(x0±ty0)/qF2]  a'±  dA  ,  (135) 

J  AB 

where  the  generalized  amplitude  function  a'±  is  given  by 

a±  =  u±af  -  ip2q  [iq(u±)2  af  +  P2y6  x\  +  0(F4)  .  (136) 


The  low-Froude-number  asymptotic  expansion  (134)-(136)  shows  that  the 
contribution  of  the  nonlinear  term  x  in  the  free-surface  boundary  condition  in  the 
generalized  amplitude  function  a'±  is  0(F2),  which  is  asymptotically  negligible  in 
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comparison  with  the  0(1)  term  u±af  given  by  the  linear  Neumann-Kelvin  theory 

except  if  u±af  =  0.  More  precisely,  Eq.  (136)  shows  that  the  low-Froude-number 

asymptotic  expansion  obtained  in  the  foregoing  for  the  wave-amplitude  function  K(t) 

immediately  yields  a  corresponding  expansion  for  the  generalized  wave-amplitude 

function  K'(t).  Specifically,  the  asymptotic  expansion  for  K(t)  given  by  Eqs.  (107)-(119) 

is  valid  if  K(t)  on  the  left  side  of  Eq.  (107)  is  replaced  by  K'(t)  and  if  the  term  * 

iq(u±)2af  in  expressions  (56)  and  (86)  for  the  second-order  amplitude  functions  Af  is 

replaced  by  the  term  iq(u±)2a2±  +  p2yox(40- 

This  generalized  asymptotic  expansion  for  the  function  K'(t)  shows  that  we  have 
k(t)  =  0(F3),  as  may  be  seen  from  Eqs.  (107),  (111)  and  (112).  This  contribution  of 
the  nonlinear  terms  in  the  free-surface  boundary  condition  must  be  compared  with 
the  function  K(t)  corresponding  to  the  linear  Neumann-Kelvin  theory.  Specifically,  we 
have  K(t)  =  O(F)  for  a  ship  form  having  a  region  of  flare,  K(t)  =  0(F2)  for  a  ship 
form  that  is  wall  sided  everywhere  and  has  either  a  bow  or  a  stern  (or  both)  that  is 
neither  cusped  nor  round,  and  K(t)  =  0(F3)  for  a  wall-sided  ship  form  with  both 
bow  and  stern  that  are  either  cusped  or  round.  The  contribution  of  the  nonlinear 
terms  in  the  free-surface  boundary  condition  is  then  asymptotically  negligible  for  ship 
forms  having  a  region  of  flare  or  a  sharp  wedge-like  bow  or/and  stern,  which  is  the 
case  for  a  majority  of  ships,  notably  naval  ships. 
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